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Intensity interferometry in noncentral heavy ion collisions 
provides access to novel information on the geometry of the 
effective pion-emitting source. We demonstrate analytically 
that, even for vanishing pair momentum, the cross terms 
and R 2 gl of the HBT correlation function in general show a 
strong first harmonic in their azimuthal dependence. The 
strength of this oscillation characterizes the tilt of the major 
axis of the spatial emission ellipsoid away from the direction 
of the beam. Event generator studies indicate that this tilt 
can be large (> 20°) at AGS energies which makes it by far 
the most significant azimuthally sensitive HBT signal at these 
energies. Moreover, transport models suggest that for pions 
this spatial tilt is directed opposite to the tilt of the directed 
flow ellipsoid in momentum space. A measurement of the az- 
imuthal dependence of the HBT cross terms R 2 t and R 2 sl thus 
probes directly the physical origin of directed pion flow. 

PACS numbers: 25.75.+r, 07.60. ly, 52.60.+h 

Two-particle momentum correlations between identi- 
cal particles are commonly used to extract space-time 
and dynamical information about the particle emitting 
source in heavy ion collisions. The basis for this inten- 
sity interferometric method is the equation jl] |J| 



C(q,K) = l 



|J d 4 xS(x,K) e iq - x \ 
|/ d*xS(x,K)\ 2 



(1) 



q = pi — P2, K = h(px +Ps), which relates the phase- 
space density S(x,K) of the source to the measured 
2-particle correlation function C(q, K). Experimental 
measurements of C are usually parametrized in terms 
of the intercept A(K) and the HBT radii iJf.-(K) by 



C(q, K) = 1 + A(K) exp 



E 

i,j=o,s,l 



q iqj RUK) . (2) 



In this Cartesian os^-system the relative momentum is 
decomposed into components parallel to the beam (I = 
longitudinal) , parallel to the transverse component of K 
(o = out), and in the remaining third direction (s = side). 

Over the last decade the experimental frontier in 
studying identical two-particle correlations was defined 
by more and more differential measurements of the HBT 
radii R%{K). At both AGS §,§ and SPS §] energies, the 
longitudinal [Kl) and transverse {K±) pair momentum 
dependence is now well-studied for central and reaction- 
plane averaged non-central collisions. Most importantly, 
these studies have led to a detailed characterization of 



the longitudinal expansion and the transverse radial flow 
of the reaction zone. The next challenge is a similarly de- 
tailed study of the i£?- (K) as a function of the azimuthal 
orientation $ of the transverse pair momentum Kj_ with 
respect to the impact parameter b in non-central colli- 
sions. This ^-dependence reveals qualitatively new infor- 
mation about the space-time structure of the source and 
yields new insights on the underlying nature of flow. In 
fact, under reasonable assumptions all 10 components of 
the source's spatial correlation tensor can be recovered. 

An azimuthally sensitive HBT analysis involves all six 
parameters R 2 j , all of which are functions of all 3 com- 
ponents of the pair momentum K± , Y and $ Q . They 
provide information about the source in space-time ac- 
cording to the following relations j7j,|| : 

R 2 {K ± , $, Y) = S n sin 2 $ + S 22 cos 2 $ - S 12 sin 2$ , 
R 2 {Kj_, $, Y) = S n cos 2 $ + S 22 sin 2 $ + S 12 sin 2$ 
-2/3 ± S Q1 cos $ - 2p ± S 02 sin $ + 0^S W , 
R 2 os {K ± , $, Y) = S 12 cos 2$ + \ {S 22 - Sii) sin 2$ 

+0xSoi sin <f> - (3±S Q2 cos $ , 
Rf(K ± , $, Y) = S 33 - 2ftS 03 + PfS 00 , 
R 2 ol (K ± ,1>,Y) = (S 13 -ASoi)cos$-/?xS, 



+ (S 23 - (3iS 02 ) sin $ + pipxSoo , 
R 2 sl (K ±: <P,Y) = (S 2 3-ASo2)cos$ 
-(S 13 -0iS O i)sm^. 



(3) 



The pair velocity (3 = K /K° arises from the on-shell con- 
straint |ij q° = q/3; it mixes spatial and temporal infor- 
mation, denotes the spatial correlation tensor 



S [us (XfiXu) j 



■Xu, (^ = 0,1,2,3) (4) 



which measures the Gaussian width in space-time of the 
emission function S(x, K) around the point of highest 
emissivity 5 M = [||: 



(x^Xv)(K) 



J d 4 x x^x v S{x, K) 



Jd 4 xS(x,K) ' 
It is the inverse of the curvature tensor B^ v 



S(x, K) w N(K) S(x, K) exp 



-x^B x v 



(5) 



(6) 



This approximation neglects non-Gaussian components 
of the emission function whose influence on the HBT 
radii can in most practical cases be neglected ||. We 
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emphasize that in (|4|) S^ v is defined in terms of Carte- 
sian coordinates in an impact parameter fixed system, in 
which x\ = x is parallel to the impact parameter b and 
X3 = z lies in the beam direction. 

The general relations (||) separate the explicit $- 
dependence of the HBT-radii (which is a consequence 
of the azimuthal rotation of the osZ-system relative 
to (x\, X2, X3)) from the implicit ^-dependence of the 
space-time widths (x^x u )(K^_, Y, $) (which reflects a 
dependent change of the shape of the effective emission 
region) |7fl. Existing studies of (||) focussed on the de- 
tailed interplay between explicit and implicit $-depen- 
dences in the HBT radii R 2 , R 2 a and R 2 os §©0. Here 
we show, however, that some of the most striking features 
are found in analyzing the ^-dependences of R 2 ol and R 2 sl 
which so far received less attention (h]] . 

The following discussion is simplified significantly by 
the important observation that the implicit $-depen- 
dence of is weak. It can be neglected relative to 
the explicit one given in (||) as long as the ^-dependence 
of space-momentum correlations in the source is small 
compared to the thermal smearing, and for K± — >■ it 
vanishes completely. Studies with the RQMD model JlJ] 
indicate that the first condition works well at least up to 
p T = 300 MeV/c for Au+Au collisions at 2 AGeV @. 
Beyond such model studies, a simple scale argument illus- 
trates why neglecting the implicit ^-dependence relative 
to the explicit one has a much wider kinematical region of 
validity than neglecting the implicit i4fj_-dependence rel- 
ative to the explicit one in (||): the latter is suppressed 
near K± =0 ((3±_ = 0), and it multiplies only space-time 
variances involving t which are numerically small in prac- 
tice. In contrast, the explicit ^-dependence in (|3|) leads 
to prefactors cos(n5>), sin(n$) oscillating between 1 and 
— 1 even for K± = , and it multiplies the numerically 
large components of S^ v . The assumption of vanishing 
implicit ^-dependence can be checked experimentally , 
and deviations can be quantified in a full harmonic anal- 
ysis given elsewhere [0,^3). Also, while it requires weak 
transverse flow, there is no such restriction on the longi- 
tudinal flow. Qualitatively, the main findings presented 
here do not depend on this assumption, but it simplifies 
our presentation and allows for a particularly intuitive 
geometric picture of the new effect discussed here. 

With this proviso, the components <S^ t „ in (|^) become 
^-independent constants which describe the same source 
being viewed from all angles We turn briefly to sym- 
metry considerations at midrapidity. Considering col- 
lisions between equal mass nuclei, it can be rigorously 
shown |l3| ] that, as a consequence of point reflection sym- 
metry around the spatial origin and mirror symmetry 
with respect to the reaction plane, five of the off-diagonal 
components (all except 5*13) oscillate symmetrically 
around zero. Coupled with the condition of vanishing 
implicit ^-dependence this implies 



S01 = , S02 = , S* 03 = , S12 = , S 23 = . (7) 

These equations and the fact that around midrapidity 
the average (3i is zero (although average (3f ^ 0) allow 
us to write the HBT radius parameters in terms of 5 
non-vanishing components only: 



R: 



\ (S u + S22) + US22 - S u ) cos 2$ , 



R 2 Q = \ (Su + S 22 ) - I (S22 - S u ) cos 2$ + f3lS o , 
R 2 os = US 22 -Su)sm2<Z> 



R 2 = S33 + f3?Soo , 
R 2 ol = S13 cos $ , 
R 2 sl = —S13 sin $ . 



(8) 

Since the ^-dependences of R 2 and R 2 explicitly sepa- 
rate S22 from Su, the emission duration Sqq can now be 
determined without additional model assumptions , 
contrary to the case of collisions at b = 0. 

Given the measured weak F-dependence of the HBT- 
radii j|"U, Eqs. (||) can be used in practice also for event 
samples which are averaged over large K-windows sym- 
metric around Y = 0. According to (||), the HBT radius 
parameters R 2 , R 2 and R 2 S all show second harmonic os- 
cillations of the same strength \ (S11— S^)- This is the 

i?o 2 2 — _ ^s/ = — R-os 2 2 rmc f° r second harmonic coef- 
ficients Q; leading deviations from this rule have been 
quantified [7|JT3]| and provide a consistency check on the 
assumption of negligible implicit ^-dependence. More 
strikingly, R 2 ol and R 2 sl display purely first harmonic os- 
cillation at midrapidity which are easier to measure. The 
expected identical amplitudes for these oscillations pro- 
vide a further consistency check on our assumptions. 

The required $-binning puts severe demands on the 
pair statistics. A first observation of azimuthally oscillat- 
ing HBT radii was reported in |l4j ]. Such measurements 
require a reasonably accurate determination of the reac- 
tion plane; a typical uncertainty of 30° reduces the first 
and second harmonics in (||) by ~ 15% and ~ 45%, re- 
spectively, but these losses can be corrected for fiplf . 

While the amplitude of the oscillations of R^R 2 , 
and R 2 S are given by the difference between the trans- 
verse source sizes in and perpendicular to the reaction 
plane, that of the oscillations of R 2 sl and R 2 ol is given by 
S\3 = (xz). Parameterizing the source by an ellipsoid, a 
nonzero S13 corresponds to a tilt in the reaction plane of 
the longitudinal major axis of the ellipsoid away from the 
beam direction. It can be characterized by a tilt angle 
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2 



S33 — Su 



Rotating the spatial correlation tensor S^ v by 9 S yields 
a purely diagonal tensor S' = Rl(0 s ) • S • R y (9 s ) whose 
eigenvalues are the squared lengths of the 3 major axes. 

We illustrate the role of the tilt angle (0) with a tilted 
Gaussian toy distribution with no space-momentum cor- 
relations: 
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S(x,K) = e- B / T exp 



x cos — z sin . 



„/2 



.72 



i 2 



a;'" y z" 
z' = a; sin 6 + z cos 6. (10) 



To avoid relativistic complications, the "temperature" T 
is kept small (20 MeV) in the following. Fig. |l| shows the 
projection of this source onto the reaction (xz) plane. 




FIG. 1. The tilt of the spatial distribution of pion emis- 
sion points projected onto the reaction (xz) plane. Contours 
show on a logarithmic scale the spatial distribution for the toy 
model ( fu5| ) with crt = 5fm/c, cr x = 4 fin, a y — 5im, cr z =7fm, 
and 6 = 25°. The shaded region is the distribution of emis- 
sion points of pions with \p z \ < 40 MeV/c in the toy model 
including longitudinal flow. See text for details. 

Using the model JIc| ) with the parameters of Fig. [I] to 
randomly generate a set of phase-space points, we con- 
structed a three-dimensional correlation function (with 
« 4x 10 5 pairs with q < 100 MeV/c) for each of eight 45°- 
wide <f> bins, according to the prescription and code of 
Pratt [fL6[ . Fitting each with the Gaussian parametriza- 
tion (|2j) yields the HBT radii presented in Fig. |2|. Treat- 
ing the ten components of S^ v as parameters, we perform 
a global fit with Eqs. (^) on these ^-dependent radii. The 
fit results are indicated by solid lines in Fig. ^. Applica- 
tion of (|) to the fit results yields 9 S = 24.6° ± 0.6°, in 
good agreement with the input value 9 = 25°. The diago- 
nal elements S'^ reproduce, within statistical errors, the 
input values for the squares of the homogeneity lengths: 
a x = 4.05±0.05 fm, a y = 4.97±0.04 fm, a z = 6.97±0.04 
fm, and o~ t = 4.78 ± 0.45 fm. (The larger uncertainty in 
at arises from the low /3j_ of the pions in our example.) 

While one may escape the effects of transverse flow 
(which may generate a ^-dependent effective source) by 
selecting pion pairs at low K±, longitudinal flow, which 
generates z — p z correlations, is generally stronger and 
cannot be cut away. Fortunately, since they are essen- 
tially orthogonal to the azimuthal dependences we are 
discussing, such correlations do not drastically alter the 
intuitive geometric picture we have discussed - the same 
source is still viewed from all angles <£>. 

As an example we added a boost-invariant longitudinal 
flow component in z-direction to our toy source (scaled 
so that the collective flow velocity at z = ±c z is equal to 
the thermal velocity), leaving the geometry unchanged. 



This results in (i) an increase in the tilt angle 6 S from 25° 
to 33° and (ii) a reduction in S' 33 from 49 fm 2 to 31 fm 2 . 
The other components S'^ vary negligibly from the sce- 
nario without flow. Familiar from the case of azimuthally 
symmetric HBT, effect (ii) is understood in terms of a re- 
duction in the length of homogeneity due to the flow M : 
HBT correlations arise from particle pairs with close-by 
momenta; the space-momentum correlations induced by 
longitudinal flow then imply that they will be close-by 
in coordinate space as well. The increased tilt is simi- 
larly understood, by examining the spatial distribution 
of emission points for pions with low p z . The shaded re- 
gion in Fig. |l| shows the effective source for pions with 
\p z \ < 40 MeV/c; it is clearly less prolate and more tilted 
than in the case of no flow (contour lines). 
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FIG. 2. Azimuthal dependence of the HBT radii from fits 
to correlation functions generated from the toy source ([To|). 
Solid lines represent a global fit to -R($) with (|3n. The value 
from 4> = 0° is replotted with an open symbol at $ = 360°. 

^-sensitive HBT studies and the measurement of spa- 
tial tilt connect for the first time the physics of directed 
flow with the space-time structure of the source JToJ . We 
consider this an essential step towards a full understand- 
ing of the phase-space dynamics of heavy-ion collisions. 
Although a detailed discussion must await a longer pa- 
per, we here shortly touch on the main physics points, 
using results from a realistic transport model. 

We performed simulations of semiperipheral Au+Au 
collisions at 2 AGeV with the RQMD (v2.3) model 
The top panel of Fig. |^ shows (p x ) - the average pion 
momentum in the reaction plane - as a function of mo- 
mentum p z along the beam axis. Qualitatively consistent 
with experimental observations |]l7f , a very weak negative 
directed flow ( "anti-flow" ) signal is observed - the aver- 
age emission ellipsoid in momentum space is tilted to a 
negative angle with respect to the beam (the direction of 
directed proton flow defines the positive direction). The 
magnitude of the collective motion (~ lOMeV/c) is small 
compared to the typical p? scale (~ 200MeV/c); hence 
thermal smearing dominates. The bottom panel shows 
that, while the spatial distribution displays a richer struc- 
ture than our toy model, it is nevertheless always charac- 
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terized by a significant positive tilt - opposite the average 
tilt in momentum space. A full correlation function ana- 
lysis of the RQMD events jl3| yields qualitatively similar 
results as those shown in Fig. 0. 




FIG. 3. RQMD simulation of pions from 2 A GeV Au+Au 
collisions at 6=3-7 fm. The top panel shows a weak p x —p z 
"anti-flow" correlation. In the bottom panel, contours of the 
spatial distribution of emission points projected onto the reac- 
tion plane show a strong tilt in the opposite direction from the 
tilt in momentum space. Superimposed arrows represent the 
average pion momentum at different values of z. Note that 
the momentum scale in z-direction is compressed for clarity. 

This bears directly on the physical causes of directed 
pion flow at these energies. Detailed transport model 
studies 1 18 1 have shown that pion reflection from (not 
absorption by) the nucleonic matter is at the root of 
directed pion flow at these energies. Focussing on the 
forward hemisphere, if absorption processes (jtNN 
AN — > NN) were dominant in producing pion flow, we 
would expect an absence of 7r emission points in the +x 
quadrant, i.e. a negative tilt in coordinate space and in 
momentum space. Since it is the point of last scattering 
(as opposed to the original point of creation) which is rel- 
evant for HBT correlations ||, it is clear that reflection 
(ttN — > A — > ttN) from flowing participant or spectator 
baryons leads to a positive tilt in coordinate space as seen 
in Fig. [3| the reflected pions "illuminate" the coordinate- 
space anisotropies of the nucleonic matter. In this simple 
picture, then, the sign of 9 S immediately distinguishes 
between these two possibilities. 

The arrows in Fig. ^| represent the average momenta 
of pions for different values of z. The resulting structure 
further underscores the importance of pion rescattering: 
Clearly, the more numerous pions from the high-density 



region around z — dominate, generating the anti-flow 
signal seen in experiment. However, pions from the more 
dilute large- 1 z | region have less opportunity for rescatter- 
ing and so retain the positive p x —p z correlation of their 
(flowing) parent A's. Similar considerations generate a 
sign change in the pion flow as the impact parameter is 
varied in transport models Jl8[ . 

In summary, for non-central collisions all ten compo- 
nents of the spatial correlation tensor are accessible 
by ^-dependent HBT measurements. Based on symme- 
try and scale considerations we argue that for low K± 
the explicit ^-dependence of Eqs. (|^) dominates. Con- 
sistency relations allow to check whether this is true in 
practice. The spatial correlation tensor S^u can then be 
extracted completely from a global fit to the six <I>-de- 
pendent HBT radii. At midrapidity, the five nonvanish- 
ing components of correspond to the four spacetime 
lengths of homogeneity and a tilt of the source in the 
reaction plane, away from the beam direction. This tilt, 
which may be quite large at AGS energies, causes strik- 
ing and relatively easily measurable first-order harmonic 
oscillations in R a i and R s i and can give a direct experi- 
mental handle on the origin of pion flow at these energies. 
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